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ABSTRACT 

The  paper  deals  with  life  distributioms  for  coheremt 
systems  of  compoaeats.  Two  major  questioae  are  discussed: 

(i)  eatimatioa  of  system  life  from  data,  or  kaowladge,  om 
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The  asymptotic  models  are  discussed  in  the  light  of 
hazard  rate  properties  of  the  limiting  distributions 
of  the  models*  A  parametric  family  of  distributions  is 
proposed  for  approximating  life  distributions  whose 
hazard  rate  is  bath-tub  shaped ,  thus  representing  a 
burn-in  period,  an  accidental  failure  period  and  a 
wear-out  period. 
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1.  Coherent  systems:  definitions  mnd  basic 

properties 


Let  s  system  consist  of  a  components.  We  shall  refer 
to  n  ns  the  size  of  the  system.  Our  interest  both  in  each 
component  and  in  the  system  is  whether  it  functions  or  it 
failed  to  function.  By  measuring  time  from  a  fixed  origin, 
we  introduce  to  the  j-th  component  the  random  life  length 
Xj  t  the  first  point  in  time  when  component  j  fails  to 
function.  At  eanh  fixed  point  x  of  time,  we  define 


[j(i> -{: 


if  Xj  >  x, 
if 


that  is,  Ij(x),  the  indicator  of  the  event  £Xj>xJ,  is 
one  if  the  j-th  component  functions  at  time  x,  and  zero 
if  it  failed  bffore  or  at  x.  When  it  does  not  lead  to 
confusion,  we  sometimes  suppress  the  variable  x  ,  and 
we  refer  to  I ^  as  the  indicator  of  the  j-th  component's 
functioning  or  failure. 

We  assume  that,  at  any  given  time  x,  the  function¬ 
ing  and  failure  states  of  the  components  uniquely  determine 
whether  the  system  functions  or  failed  to  function.  In 
other  words,  there  is  a  function  S  of  n  variables, 
which  takes  the  values  1  and  0  only  and  such  that 
S  ■  Sx  -  Sd^x),  I2(x),...,  Im(x))  equals  one  if  the 
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system  functions  at  time  x  and  zero  if  it  felled  before 
or  at  x.  We  cell  S  the  structure  ( fuactloe)  of  the 
system.  The  time  T  to  the  first  failure  of  the  system 
is  celled  the  system  life,  or  oa  some  occasions,  the 
failure  time. 

Notice  that  each  variable  of  S  ie  either  zero  or 
one.  Whea  a  variable  of  S  is  iaoreased  from  zero  to 
ome  (that  isf  a  failed  compoaeat  is  replaced  by  a  func¬ 
tioning  oae),  we  assume  that  it  could  mot  have  a  negative 
effect  oa  S.  Ia  other  words,  if  S  *  1  whea  I^(x)  =  0 
(that  is,  the  system  functions  although  the  j-th  compoaeat 
failed),  then  S  =  1  also  when  Ij(x)  ■  1  (that  is,  the 
system  also  fumctloas  when  the  j-th  component  functions), 
assuming  that  there  is  no  change  in  the  states  of  the 
other  components*  When  this  evident  requirement  is  satis- 
fled,  we  call  the  system  (or  the  structure  S)  coherent* 

The  mathematical  equivalent  of  the  above  described  property 
is  that  the  structure  function  S  is  momdecreasing  ia 
each  of  its  variables*  We  shall  also  assume  that  eaoh 
component  serves  some  purpose  in  the  sense  that  S  is 
mot  constant  in  any  of  its  variables* 

There  are  two  special  systems  whioh  will  play  specific 
roles  in  the  sequel  and  which  we  define  below* 

Series  system.  If  a  system  fails  as  soon  as  one  of 
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its  components  fails,  we  call  it  a  series  system.  For 
such  a  system,  S  =  1  if,  aad  omly  if,  each  1^  *  1. 

Heace 

S(I1»  In)  =  miaUj,  I2f-,  In)» 

Similarly,  by  defiaitioa, 

T  —  min(X^t  Xg,***,  X^) • 

Parallel  system.  We  call  a  system  parallel,  if  it 
functioas  as  loag  as  oae  of  its  components  fuactioas.  It 
thus  follows  that 

S(Ilt  IB)  *  maT( I-^ ,  » •  •  • »  1^) 

aad 

T  3  nax(X^,  1^)* 

Siaoe  each  Ij  is  either  0  or  1, 

•ia(Ij,  Ig»  •  •  •  *  Ig)  "  ^1^2*  *  *  *a  *  CO 

aad 

1  -  max(Ilf  Im)  -  (1  -  IjHl  -  I2)...(l  -  IB),  (*■) 

where,  oa  the  right  hand  sides,  ordlaary  multiplioatioa 
Is  aoplled.  This  flexibility,  of  oourse,  is  aot  valid  for 
any  other  variables  (other  than  those  taking  the  values 
0  aad  1  only),  aad  thus  an  expression  for  S  may  aot 
he  a  qulek  tool  for  finding  T  as  a  function  of  the  X^  • 


\ 

) 


Tke  identities  (1)  and  (2)  also  indicate  tkat  one  can 
kave  several  expressions  for  S  (wkich  are  naturally 
Identical  In  value),  and  tke  same  Is  true  tot  T.  Tkere 
is,  kowever,  one  expression  for  S,  In  wkich  only  tke 
functions  max  and  min  occur,  and  in  this  form,  if 
is  replaced  by  X^,  tken  S  transforms  into  T. 

For  formulating  tkis  special  representation,  let  us 
start  with  an  example. 

k  -  out  -»  of  -  n  system.  A  system  is  called  a 
k-out-of-n  system,  if  its  size  is  n  and  if  it  func¬ 
tions  if,  and  only  if,  at  least  k  of  its  components 
function  (1  k  n).  Wltk  tke  delta  function 


A(l,i)  -  j* 


1  f 

otherwise. 


evidently 

S(l^i  I2***'*  ^  (^1  ^2  *•••*’  ^n* 

Tke  definition  also  yields  tkat 


n-k+lsn  * 


wkere  Xr:m  is  tke  r-th  order  statistic  in  tke  sequenoe 


Xj,  Xg , . . . ,  x^  (Xisa<  X2  ^ 


£•  Xmsm  )•  low,  einoe 


L 


#  •  e 
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X*-k+lsm  =  xl2 . 

A(W"  .+In,k)  =  max £mia( I ^  ,  I±  ,...,  IA  )J  ^ 

12  k 


where  l£i1<l2  c  . and  the  maximum  ie  takem 
over  all  the  possible  choices  of  the  k  subscripts  i ^ , 

1 .5 j  5k,  we  expressed  S  aad  T  by  the  same  fuactioa  of 
the  variables  I ^  aad  ,  15j<n,  respectively. 

The  expressioms  for  S  aad  T  as  the  maximum  of 
some  variables  shows  that  a  k-out— of-a  system  caa  also 
be  viewed  as  a  parallel  system.  Iadeed,  if  we  fix  the  in¬ 
tegers  1 5  i.^  ig  <  •  •  •  <  ifc  5  a  aad  we  coaneot  these  k 
components  into  a  series  system,  then  the  parallel  system 
whose  components  are  the  just  constructed  j  series 
systems  is  equivalent  to  the  original  k-out-of-n  system 
la  that  their  structure  functions  aad  their  system  lifes 
are  identical. 

The  fact  that,  by  suitable  grouping  of  components,  a 
system  can  be  reduced  to  a  parallel  system  is  not  speolal 
to  the  k-out— of-n  system.  Iadeed,  an  arbitrary  coherent 
system  oaa  be  reduoed  to  a  parallel  (or  a  series)  system. 
For  proving  this  fact,  we  introduce  the  following  concepts. 

Minimal  path  sets.  A  path  set  of  a  coherent  system  is 
a  set  C  of  components  such  that  if  each  member  of  C 
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functions  then  the  system  functions.  A  path  set  Is  mini¬ 
mal,  if  the  removal  of  a  single  element  from  it  results  in 
its  failure  to  remain  a  path  set. 

Minimal  cut  sets.  A  set  C  of  components  is  called 
a  cut  set  if  the  system  fails  whenever  each  member  of  C 
fails.  A  cut  set  is  minimal  if  no  element  can  be  removed 
from  it  without  violating  its  cut  set  property. 

Now,  a  coherent  system  which  is  capable  of  functioning 
is  necessarily  boti,a  pa$h  set  and  a  cut  set.  Therefore, 
every  such  coherent  system  has  at  least  one  minimal  path 
set  and  one  minimal  cut  set.  Let  A^,  A^,...,  Ap  and  C^, 

C 2»...»  Cm  be  the  distinct  minimal  path  and  cut  sets, 
respectively,  for  a  given  coherent  system.  Them^by  defini¬ 
tion, 

S(I, ,  10 .  I_)  »  max  5  min  Ii  ( 

1  *  m  I5t<p  L  3€At  3  J 

*  min  4  max  L  r  , 

i<t<m  1 3 ect  3  J 


and 


T  = 


max  A  min 
l^t  <p  (j£At 


min  s  max  X.f. 
l<i<n  «3€Ct  33 


Putting 


U 


t 


min  X.  , 
d€At  3 


3  €Ct 


(5) 


we  obtained 


T  =  wax  0.  =  aim  V.  . 
l<t£p  x  Ut<a  x 

la  other  words,  the  eyetew  life  of  an  arbitrary  coherent 
system  Is  an  extreme  value  (maximum  or  minimum)  of  a 
suitably  chosen  sequenoe  of  random  variables.  Its  distrib¬ 
ution  function,  which  we  call  a  failure  time  distribution 
(for  coherent  systems),  is  therefore  an  extreme  value 
distribution  In  some  appropriate  model.  It  should  be 
emphasised  that  independence  and  stationarity  cannot  be 
assumed  here  even  if  the  orogonal  components  are  believed 
to  function  Independently  (such  a  case  is  rare  though, 
see  the  last  hut  one  paragraph  in  the  next  section  in 
this  regard).  Because  the  minimal  path  sets  (or  the 
minimal  cut  sets  C^)  are  not  disjoint,  the  random  variables 
and  of  (3)  are  strongly  dependent.  Their  dependence 

is  determined  by  the  structure  of  the  underlying  system 
and  by  the  dependence  of  the  original  componants.  Hence, 
their  dependence  is  never  a  matter  of  arbitrary  assumptions. 
Thus  the  study  of  these  distributions  is  an  integral  part 
of  the  theory  of  extremes  for  dependent  models.  While 
such  a  theory  is  well  developed  in  Chapter  3  of  the  present 
author's  booh,  Galambos  (1978),  we  shall  discuss  some  points 
of  this  theory  as  they  relate  to  failure  time  distribu¬ 
tions.  Although  in  the  present  paragraph  we  defined  failure 
time  distributions  in  terms  of  coherent  systems,  much  of 
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what  le  to  be  said  aleo  applies  to  failure  times  of  dans, 
when  failure  ie  caused  by  high  floode;  of  air  auality, 
when  failure  is  defiaed  by  the  fact  that  some  specific 
pollutant  exceeds  a  given  level  of  concentration,  and  to 
other  diverse  fields.  We  shall,  however,  remain  in  our 
discussion  at  coherent  systems,  and  we  use  the  notations 
of  the  present  section  throughout  the  paper. 

The  concepts  and  relations  discussed  here  first 
appeared  in  the  paper  Bimb&um  et  al.  (1961).  Later 
extensions  of  these  relations  are  nicely  brought  together 
into  a  theory  by  Barlow  and  Proschan  (1975).  Our  discussion 
here  does  not  overlap  with  the  content  of  this  book. 

2.  Asymptotic  extreme  value  distributions 
as  failure  time  distributions 

We  have  seen  that  T  can  always  be  expressed  as  the 
maximum  or  the  minimum  of  some  random  variables.  Since  in 
these  representations  either  p  or  m  is  large  for  a  system 
with  a  large  number  n  of  components  (recall  that  we 
agreed  not  to  ooneider  nonessential  components  defined  as 
components  in  whose  indicators  the  structure  function  S 
is  constant),  we  assume  for  the  general  discussion  that 

T  ■  max  U 
1  ^  t  <  p 

with  p  large.  Now,  if  the  distributions  and  the  inter- 
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dependence  of  the  are  known  them  the  distribution  of 
T  is  uniquely  determined.  Hence,  only  that  cnee  presents 
a  problem  when  either  the  distribution  functions  F^x) 

■  P(U^.<Cx),  1  it  £p,  or  the  dependence  relation  of  the 
are  unknown.  In  such  cases,  an  approximation  becomes 
necessary.  It  is  shown  on  p.  90  of  Galambos  (1978)  (and 
further  explored  in  Galambos  (1981))  that  a  reliable 
approximation  to  the  distribution  of  T  cannot  be  obtained 
through  an  approximation  to  F^(x).  Rather,  one  should 
develop  a  dependent  model  for  the  ,  evaluate  the 
possible  limiting  distributions  for  the  maximum  in  that 
model,  and  one  of  these  possibilities  is  to  be  applied 
as  an  approximation  to  the  distribution  of  T. 

There  are  a  number  of  dependent  models  for  which 
the  mathematical  results  are  at  an  advanced  stage  (although 
far  from  complete).  These  are  described  in  Chapter  3  of 
the  mentioned  book  of  the  present  author,  from  which  we 
yuote  the  following  results. 

(i)  Approximation  by  the  olassical  model  (the 
are  Independent  and  identically  distributed) 
is  ver  rarely  justified,  but  when  it  As  appli¬ 
cable,  then  the  failure  time  distribution  is 
Weibull. 

(li)  The  assumption  that  the  form  a  sequence  of 
exchangeable  random  variables  is  mathematically 
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always  justified*  But,  because  of  this  gener¬ 
ality,  the  possible  limit  laws  for  the  maximum 
form  a  very  large  family.  The  invest igatioa  of 
the  properties  of  some  subfamilies  of  these 
distributions  would  be  a  very  important  task* 

(it  should  be  remarked  that,  contrary  to  the 
claim  in  the  paper  Zidek,  et  al.  (1979),  the 
sequence  cannot,  in  general,  be  considered 
as  a  segment  of  an  infinite  sequence  of  ex¬ 
changeable  variables;  only  finite  exchangeabil¬ 
ity  is  justified.) 

(ill)  There  is  a  general  dependent  model,  in  which 

the  possible  limit  laws  for  the  maximum  coincide 
with  the  distributions  whose  hazard  rate  func¬ 
tion  idefined  below)  is  monotonio  (see  Sec¬ 
tions  3*9  and  3*10  in  the  quoted  book).  Beoauae 
of  the  significance  of  this  result  to  engineers, 
we  discuss  this  conclusion  in  more  detail* 

Engineers  have  recognized  for  a  long  time  the  impor¬ 
tance  of  the  hazard  rate  function  of  a  failure  time  dis¬ 
tribution.  Let  us  first  give  the  definition  of  hazard 
rate. 

Let  X£0  be  a  random  variable  with  distribution 
function  P(x)  and  with  density  function  f(x)  *  ?*(x).  Then 
the  hazard  rate  r( x)  of  X  is  defined  by  the  limit  relation 
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r( x)  *  lim  — —  P(X<Tx  +  <Tx  I  X  £x) . 
^x=0  ^  *  « 


Am  easy  calculation  yields  from  this  limit  relation 


fUJ 

r(x)  =  - -  . 

1  -  FU) 


($ ) 


While  this  latter  form  is  a  convenient  formula  for 
actually  calculating  r(x),  the  definition  of  r(x)  is 
what  makes  it  applicable.  When  X  represents  a  random 
life,  then  the  conditional  probability  P(X<  x  +  «Tx  [  Xix) 
expresses  the  probability  of  X's  failing  in  the  short 
interval  ( x,  x  +  sTx),  given  that  X  has  survived  beyond 
x.  It  is  apparent  to  an  engineer  that  a  new  system  of 
eomponents  whose  life  is  represented  by  X  may  have 
eome  positive  probability  of  failing  immediately  after 
production,  but  this  probability  decreases  as  time 
passes  (bum- in  period).  On  the  other  hand,  an  old  system 
of  components  is  more  and  more  likely  to  fail  as  time 
passes  (ageing  or  wear-out  period).  Between  the  burn-in 
and  wear-out  periods,  for  most  systems,  it  is  accepted  N 

(met  vV<.r>*t^L  -fw-ttMrc  pe.r»»w/ 

that  only  accidents  may  cause  fnilureY'Tke  stochastic 
definition  of  aocidents  is  either  by  constant  hazard  rates 
or  by  the  laok  of  memory  property 

P(X£  x  +  y  |  X£  x)  *  P(X£y). 

Sea  Section  1.5  in  Galambos  and  Kotz  (1978)  to  see  that 
these  two  seemingly  different  definitions  are  equivalent. 


14. 


Noticing  that  oomstaat  hazard  rate  caa  also  he  viewed 
ae  both  mondecr  easing  and  nominor  easing,  that  ie  mono tonic, 
we  conclude  that  the  life  length  of  a  system  can  be  split 
into  three  periods,  in  each  of  which  the  hazard  rate  of 
the  life  distribution  is  monotonic.  It  is  very  pleasing 
to  see  that  the  mathematical  theory  through  structure 
functions  and  extreme  value  theory  led  to  a  similar 
conclusion. 

Returning  to  the  simple  argument  of  the  preceding 
paragraph,  we  can  also  find  it  reasonable  that  the  smaller 
is  the  change  in  the  hazard  rate  the  closer  we  are  to  tha» 
aocldental  failure  period  (that  is,  sharper  decreases 
occur  at  the  beginning  of  the  burn-in  period  than  later 
in  this  same  period  and  sharper  and  sharper  increases  in 
the  hazard  rate  are  experienced  during  the  wear-out  period 
as  tine  passes).  ®his,  when  translated  into  mathematical 
terns,  implies  that  the  hazard  rate  function  is  convex, 
whloh,  together  with  the  three  periods  discussed  earlier, 
is  frequently  referred  to  as  the  hazard  rate  funotion  is 
bath-tub  shaped. 

A  good  approximation  to  most  empirical  results  oam 
be  obtained  by  the  following  bath  tub  shaped  hazard  rate 
funotion: 


A(x  -  a)2  +  c 

if 

0  <  x<a 

c 

if 

a  <  x^b 

B(x  -  b)2  ♦  c 

if 

x£b 

r(x; A,a,c,B,b) 
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The  coast ant  c  >  0  represents  the  constant  hazard  rate 
during  the  accidental  failure  period  (a,b),  while  the 
btirn-ln  period  lasts  upto  the  point  a  and  the  wear-out 
period  starts  at  b  .  The  constants  A>0  and  B>0 
are  shape  parameters.  The  choice  of  a  quadratic  function 
is  arbitrary  here,  but  when  an  approx imar ion  is  to  be 
found  on  an  empirical  basis  rather  than  by  theoretical 
arguments,  then  the  shape  parameters  give  sufficient 
flexibility  to  get  a  good  approximation.  We  add  here 
that  there  is  no  result  regarding  the  statistical  analysis 
of  the  distribution  whose  hazard  rate  is  r(x;A,atc,£,b) . 
Its  study  would  be  very  useful. 

There  is,  of  course,  a  well  determined  relation 
between  a  distribution  function  and  its  hazard  rate 
function.  Namely,  if  we  write  (5)  in  the  form 

d{ln(l  -  P(x))J 

r(x)  *  -  -  , 

dx 


then  we  immediately  see  that  P(x)  unifuely  determines 
r(x).  On  the  other  hand,  we  get  by  integration  from  the 


above  relation. 


f  * 

?(x)  *  1  -  expj -Jr(t>dt|,  x>0. 

o  J 

We  thus  see  that  when  the  hazard  rate  is  constant,  then 


the  distribution  function  is  exponential*  In  particular, 


during  the  accidental  failure  period,  failure  time  dlstribu— 


tioas  are  always  exponential.  This  observation  leads 
us  to  the  following  important  ooaclueloa.  Assume  mow 
that  all  components  as  well  as  the  system  achieved  the 
aocideatal  failure  period.  This  means,  that  both  the 
oompoaents  and  the  system  have  exponential  failure 
distributions  (for  a  certain  period  of  time  only).  It 
them  follows  from  a  result  of  Esary  et  al.  (1971)  that 
the  components,  with  the  exception  of  series  systems, 
are  not  stochastically  independent.  That  is,  one  cannot 
construct  a  single  structure  other  than  a  series  system 
in  which  the  components  would  function  independent: 

(and  whioh  system  would  aohieve  an  accidental  failure 
period).  This  Is  a  very  important  conclusion  because 
several  estimates  on  reliability  are  developed  in  the 
literature  under  the  assumption  that  the  components  are 
Independent. 

Finally,  we  remark  that  there  are  a  number  of 
characterisation  theorems  for  exponent iality  (see  the 
book  Galambos  sad  Eots  (1978))  which  can  be  used  for 
testing  whether  a  system  Is  in  its  accidental  failure 
period.  In  most  oases,  those  limited  characterisation 
theorems  are  sufficient  when  one  assumes  a  priori  that 
the  underlying  distribution  is  of  aonotonic  hasard  rats 
A  typical  result  of  this  nature  oan  be  found  in  Ahsa- 
nullah  (1977). 


17. 


3.  Estimates  on  failure  time 
distributioae 

We  have  emphasized  ia  the  previous  section  that 
components  for  most  structures  cannot  function  indepen¬ 
dently  of  each  other.  At  the  same  time,  we  may  know 
exactly  the  distribution  of  component  lives,  mainly 
through  characterization  theorems.  This  leads  us  to  the 
problem  of  estimating  failure  tine  distributions  by 
the  distributions  of  component  lives  under  some  assump¬ 
tion  of  dependence  of  the  components. 

We  use  the  path  set  and  out  set  decompositions, 
in  view  of  which  structure  life  is  an  extreme  of 
"component  lives"  (where  "component"  is  either  a  minimal 
path  set  or  a  minimal  cut  set).  Through  this  approach, 
the  mathematical  problem  is  the  estimation  of  the 
distribution  function  H (x)  of 

H 

T  *  max( U^,  »•••#  Up) 

under  some  form  of  dependence  of  the  and  under  the 
assumption  that  the  distribution  functions  F^(z)  • 
P(Ut<  x)  are  known. 

There  is  one  ooncept  of  dependence,  the  so  called 
association  of  random  variables,  for  whloh  there  is  an 
extensive  literature  with  reliability  emphasis  (see 
Barlow  and  Proechan  (1975)  and  Natvig  (1980)).  However, 


these  works  deal  with  estimating  E(T)  la  terns  of  E(Ut), 
l^t  5p,  rather  than  giving  estimates  on  HJix)  • 

Since  we  deal  with  distributions,  we  express  de¬ 
pendence  through  distributional  assumptions.  The  simplest 
distributional  assumption  is, of  course,  when  only  bivariate 
distributions  are  involved.  Por  simplicity,  we  Introduce 
the  notations  =  A^(x)  * 


51, p(x)  *  Sl,p  *  ^  P*Aj*» 

j  =  4  _ 

52, p(x)  *  S2,p  *  Z*-  PUjAAj), 

and  we  let  m^  *  mp(x)  stand  for  the  number  of  those  A^ 
which  oocur.  Then  H^x)  *  PCn^  *  °)»  Sl,p  *  E<V 
and  2S2  p  +  si,p  *  This  latter  meaning  of 

and  S0  makes  them  appealing  to  the  applied  statistician, 
while  their  original  definition  is  the  useful  fora  in 
mathematical  arguments.  It  is  slightly  more  convenient 
to  state  results  for  1  -  HjLx)  *  than  for  H#(x) 

Itself,  and  we  shall  do  so  below.  let  us  oonsider  esti¬ 
mates  of  the  fora 

»  Sl,p  +  *  S2,p  £  «Vl)  ^  °  S1.P  *  *  S*.P  ’  <6> 

where  a,  b,  o  and  d  are  constants  (which,  in  prinoiple, 
may  depend  on  x  suppressed  in  all  notations). 
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The  beat  lower  bound  in  (6)  la  known  (Kwerel  (1975) 
and  Galamboa  (1977)),  according  to  which  a  and  b 
ahould  be  of  the  form  a  *  2/(k+l)  and  b  ■  -2/k(k+l), 
where  l:'kv.p  in  an  integer.  It  ie  then  easy  to  find 
the  optimal  k  which  equala  2i>2  p~]  ♦  1»  where  [yj 

signifies  the  integer  pert  of  y.  (Notice  that  k  *  l  yields 
the  classical  estimate  by  the  method  of  inclusion  and 
exclusion.)  For  the  upper  bound  in  (6),  only  partial 
results  are  available.  The  best  known  result  is  c  ®  1 
and  d  «  -2/p  (Kounias  (1968)  and  Galambos  (1975)). 

Before  proceeding  with  the  discussion  of  the  estimates 
in  (6),  notice  the  following  important  fact.  The  results 
quoted  in  the  previous  paragraph  are  such  that  the  coef¬ 
ficients  a,b,c  and  d  do  not  depend  on  x.  Heaoe,  they 
remain  valid  if  we  redefine  A^.  Now,  if 
then 

K  -  °]-K<Xl'  U2<x2»‘--»  up<xp  l  * 

and  thus  (6)  provides  estimates  on  the  p— variate  distribu¬ 
tion  of  the  Dj  in  terms  of  univariate  and  bivariate  mar¬ 
ginals.  These  inequalities  should  be  taken  into  account 
when  one  is  interested  in  constructing  multivariate  dis¬ 
tributions  with  given  (univariate  and  bivariate)  marginals. 

Let  us  return  to  (6).  Since  P(m,p>l)  is  related  to 
the  distribution  of  the  maximum  of  the  U one  would  like 
to  get  suitable  extensions  of  (6)  to  P(mp£  r),  r}l, 
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which  is  relevant  for  the  distribution  of  the  (p-r+l)-st 
order  statistic  of  the  0^  .  Two  methods  of  proof  of  (6) 
lead  to  interesting  results  in  this  direction.  One  proof, 
which  roughly  says  that  (6)  is  valid  for  arbitrary  (de¬ 
pendent)  sequence  of  0^  if  it  is  valid  when  the  are 
i.l.d.,  implies  that  a  set  of  coefficients  in  (6)  deter¬ 
mines  a  set  of  coefficients  for  estimating  P(mp>r) 
in  the  form  of  (6)  (Galambos  and  Mucci  (1980)).  The 
other  proof,  introduced  in  Galambos  (1977),  provides  a 
technique  which  can  be  used  with  success  in  more  general 
cases  than  (6)  (for  example,  to  P(mp>r)  for  any  r>l, 
and  with  bounds  not  necessarily  linear).  On  this  line 
of  extension  of  (6),  we  mention  Sathe  et  al.  (1980). 

Por  earlier  results  on  inequalities  of  the  nature  of  (6), 
see  the  survey  at  the  end  of  Chapter  1  of  Galambos  (1978). 
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